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SIMULATION OF A COMBINED-CYCLE ENGINE 


Abstract 


A FORTRAN computer program has been developed to simulate the 
performance of combined— eye le engines. These engines combine features 
of both gas turbine and reciprocating engines. The computer program described 
and documented in this report can simulate both design point and off-design 
operation. Widely varying engine configurations can be evaluated for 
their power, performance and efficiency as well as the influence of altitude 
and air speed. Although the program was developed to simulate aircraft 
engines, it can be used with equal success for stationary and automotive 
applications . 
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SIMULATION OF A COMBINED-CYCLE ENGINE 


1 Introduction 

As part of a project to develop a high-efficiency engine for use in 
helicopters, the U.S. Army Aviation Systems Command authorized the 
development of a simulation capability for estimating the off-design 
performance of combined-cycle engines. These are engines that have features 
of both gas turbine and reciprocating engines. This document is the final 
report for the project to develop this simulation capability. 

During a Summer Faculty Fellowship in 1988, the author developed a 
computer program for simulating two-stroke diesel engines. This program was 
documented in reference [1]. While this program allowed prediction of 
diesel engine performance it required knowledge of the intake air supply 
pressure and temperature as well as the exhaust back pressure. This 
information can be specified at a design point operating condition but at 
an off-design condition, it is usually unknown. In order to predict 
off-design performance, it is necessary to incorporate the characteristics 
of the compressors and turbines used to supply the air to the engine. 

The Navy-NASA Engine Program (NNEPEQ) is an existing computer program that 
allows the off-design behavior of turbomachinery to be predicted and also 
solves the resulting matching problem [2,3]. The principle objective of the 
project described in this report was to combine the two-stroke diesel engine 
simulation program with NNEPEQ to obtain a combined-cycle engine simulation 
program. 

This report is suitable as a user's manual for the combined-cycle engine 
simulation program. Since the two programs that comprise the simulation are 
documented elsewhere, details of these programs will not be provided. Only 
those features of the combined program that are new, or not discussed 
elsewhere, will be documented here. 

The proposal for this project stated that in addition to the integration 
of the two simulation programs, other improvements would be made. These 
include upgrading the air scavenging model in the diesel program, 
investigating alternative compressor maps in NNEPEQ, incorporating second law 
analysis into the diesel program and allowing the use of crankshaft offset 
with opposed piston engines. While all of these features are addressed in 
this report, only the improved scavenging model has been incorporated into 
the combined cycle program. The other features are provided in special 
purpose programs. This requirement came about because of the need to 
dramatically reduce the execution time of the diesel simulation program. 

The version of the program described in reference [1] requires about 2-8 
minutes of CPU time on a VAX 11/785. Since NNEPEQ will execute this 
program approximately 100 times for each operating condition and a typical 
run of NNEPEQ involves one design point run and 8-12 off-design point 
conditions, it is obvious that CPU times had to be reduced. Program 
modifications were made that reduced the execution time by a factor of 20 with 
only minimal changes in output. The combined program is currently being run 
on a DEC 3100 workstation and requires about 3 minutes for each operating 
condition. 

This report will describe the improvements made to the diesel simulation, 
especially the air scavenging model. Then the features of the combined program 



will be presented as well as the input and output from the program. These 
features will be illustrated with two sample cases. Following the sections 
on the combined-cycle engine simulation program, the other aspects of the 
contract work statement will be discussed. An appendix is provided describing 
a computer program for doing second law calculations. Another appendix shows 
how experimental data for a small centrifugal compressor was entered into the 
program and compared with the standard maps used by the program. Finally, the 
last appendix contains a subroutine for calculating diesel engine cylinder 
volumes for opposed piston engines with crankshaft phase offset. 
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2 Structure of the Combined-Cycle Simulation, Computer Program 


Figure 2.1 shows a schematic of the combined-cycle engine simulation 
computer program. The main program, called CCE, is very similar to the main 
program in NNEPEQ described in references [2,3] . However, it contains the 
call to a subroutine INPUT1 that reads the file containing the specifications 
of the diesel engine. The block of subroutines identified as NNEPLIB includes 
all the subroutines originally associated with NNEPEQ. Only a few have been 
modified to allow the simulation of a combined-cycle engine. CCE can still be 
run as a gas turbine simulation program and retains all the features of NNEPEQ. 

The diesel portion of the simulation has been incorporated into the gas 
turbine simulation as a burner. The subroutine DBURNR from NNEPEQ was modified 
so that if a certain parameter is non-zero, the program flow branches out of 
DBURNR and calls the DIESEL subroutine. 

There is a mis-match between the requirements of the turbomachinery 
portion of CCE and the diesel program. CCE provides each of the elements 
of the gas turbine engine with an incoming pressure, temperature and mass 
flow rate. It then calculates an exit pressure that corresponds to these 
conditions. However, the diesel program expects to receive an incoming 
pressure and temperature and an exit pressure. It then calculates the mass 
flow rate for these conditions. Due to this problem it was necessary to 
develop an interface subroutine that matches the information supplied by CCE 
to that required by the diesel program. This interfacing subroutine is 
called DIESEL and it calls the diesel simulation program DSL2 with estimated 
values of exhaust back pressure. It iteratively varies this pressure until 
the mass flow rate through the diesel matches that required by CCE. 

The block of subroutines identified as ELIB includes the subroutines 
required by DSL2 . These subroutines are separate from and do not depend on 
those present in NNEPLIB. 
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Figure 2.1 Structure of the Combined-Cycle Engine (CCE) Simulation Program 
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3 Modif icat ions to the Diesel Simulation Program 


3.1 Air Scavenging Models 

The original version of the diesel simulation computer program 
described in reference [l], has a very basic scavenging model that is not 
flexible enough to accomodate differences in the scavenging characteristics of 
various engines. One of the objectives of this project was to implement a 
more sophisticated scavenging model that would provide more accurate 
predictions of air flow requirements in two-stroke engines. 

Three simple models for scavenging can be developed from the limiting 
cases corresponding to phenomena in two-stroke engines. One case, called 
"perfect scavenging," assumes the incoming air is separated from the exhaust 
gases by a membrane that prevents mixing. The incoming air pushes the exhaust 
out ahead of it and when enough air has been added to fill the cylinder, the 
exhaust gases have been completely expelled. Obviously, this is an ideal case 
and is never ' accomplished in a real engine. 

The second limiting case is "perfect mixing." In this case, the incoming 
air is assumed to mix instantaneously with the product gases in the cylinder 
as it enters the cylinder. When the exhaust process begins most of the 
expelled gases are combustion products but as the process continues, an 
increasing fraction is fresh charge. 

The final case corresponds to complete "short circuiting." This means 
that the air entering through the intake ports passes directly to the exhaust 
ports without expelling any exhaust products with it. This would be a highly 
undesirable form of scavenging since no fresh charge is retained for the 
next engine cycle. 

The original engine simulation program uses a scavenging model that 
assumes perfect mixing. While this model will accurately predict power 
and efficiency in diesel engines, it underpredicts the air flow requirements 
of the engine. Most two-stroke engines produced today have scavenging 
characteristics somewhere between perfect mixing and perfect scavenging. A 
scavenging model that takes advantage of this fact is the "mixing-displacement 
scavenging" model developed by Benson and Brandham [A]. The details of the 
model are described below. 

3.2 Mixing-Displacement Scavenging Model 

This model assumes that the cylinder is divided into two zones. One 
zone is called the mixing zone and the other the displacement zone. The ratio 
of the initial volume of the displacement zone to the total cylinder volume is 
called the displacement factor. At the point where both the intake and 
exhaust ports are open simultaneously, the start of the scavenging period, the 
exhaust gases leave the cylinder from the displacement zone following the 
commonly used rule based on isentropic flow of ideal gases through a reduced 
area section due to a pressure difference. Non-ideal behavior is incorporated 
with a discharge coefficient. Incoming air flows through the intake port 
following a similar relationship and enters the mixing zone where it is assumed 
to mix perfectly with the contents of that zone. The scavenging process 
continues until the contents of the displacement zone are completely exhausted 
and the mixing zone completely fills the cylinder. At this point, the products 
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leaving the cylinder are drawn from the mixing zone until the exhaust port 
closes. This model is a combination of a perfect mixing and perfect 
scavenging model. The displacement factor determines the relative weighting 
between the two models. This factor is chosen to provide scavenging 
characteristics that match experimental data, if available. The program 
currently uses a value of 0.55 for the displacement factor as this gives 
results typical of current production engines. 

The diesel simulation program uses a single zone combustion model. This 
model is used up to the point where the intake port opens and then the single 
zone is split into the two zones described above. Under almost all conditions, 
the model reverts to a single zone model during the scavenging process because 
the displacement zone is entirely exhausted. It is conceivable that under 
very high speed conditions or when the pressure drop across the engine is very 
low, the displacement zone might not be entirely exhausted before the exhaust 
port closes. In this case, the two zones are instantaneously mixed at the 
point of exhaust port closure and the calculation continues with the single 
zone model. 

Use of the mixing-displacement scavenging model requires no input from 
the user. It is automatically used when the program is used to simulate a 
two-stroke engine. 

3.3 Change in Gas Property Routine 

The diesel simulation routine described in reference [l] assumes that the 
contents of the engine cylinder are equilibrium combustion products at all 
times. The computer program used a subroutine, PERX, that was based on the 
subroutine PER developed by Olikara and Borman [5]. However, in a time study 
of the diesel simulation it was determined that the program was spending about 
807. of its time calculating properties. To reduce the amount of time spent in 
property evaluation a new property routine was developed that assumes the 
products consist of oxygen, nitrogen, carbon dioxide and water vapor 
under lean conditions and nitrogen, carbon dioxide, carbon monoxide, hydrogen 
and water vapor under rich conditions. This change in property routine 
resulted in a reduction in CPU time from 465 seconds to 86 seconds on a 
VAX 11/785. 

3.4 Change in Differential Equation Integrator 

The final modification to the program was to change the differential 
equation integrator used to solve the time-based equations for pressure, 
temperature, etc. in the cylinder. The original version of the program used 
DVERK , a subroutine from the commercial International Mathematical Subroutine 
Library (IMSL). This program uses the Runga-Kutta numerical procedure which 
is highly accurate but not very efficient. This subroutine also has the 
disadvantage that the program can only be run on a computer that is licensed 
to use the IMSL programs. To eliminate these problems, the subroutine LSODE 
was substituted for DVERK. This subroutine is publicly available so it can 
be easily transported to other computers. In addition, it uses an implicit 
numerical procedure that is very efficient [6]. CPU time was reduced from 86 s 
11 seconds using this routine. Accuracy was such, that no significant 
difference could be detected between results obtained with this routine and 
those derived using the IMSL program. 
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4 Inputs to the Combined Cycle Engine Simulation Program 

Input to the combined-cycle engine simulation program is provided by two 
data sets. The first data set corresponds very closely to the NAMELIST data 
input file used by NNEPEQ and completely characterizes the turbomachinery 
associated with the engine. The only modifications to this file are that if a 
diesel engine is included in the engine design, the engine fueling rate and 
speed are input through this data set. A conventional gas turbine engine data 
set can be provided to the program and it will run normally and provide the 
same results as NNEPEQ. 

The second input data set provides the specifications of the diesel 
engine. This data set is a slightly modified version of the data set used by 
the original two-stroke diesel engine program [ll. 

The turbomachinery data set is read from a file named "TURBO. INP. The 

diesel data set is read from a file named "DSL. INP". A third data file 
contains the turbomachinery maps. A default set of maps is provided with 
NNEPEQ and they were used for the example cases described later. The map file 
is read from logical unit 12. This means that an ASSIGNMENT statement must be 
used to identify the actual data file name as unit 12. 

The turbomachinery input data file is relatively complex and provides a 
great deal of flexibility in specifying engine configuration and component 
specifications. Since it is essentially identical to the input data file used 
by NNEPEQ, and is well-documented elsewhere [2,3], only the general features 
of the input file will be presented here. This report assumes the reader has 
a working knowledge of the NNEPEQ program. 

The first line in the turbomachinery NAMELIST is a title for the cases 
and can be up to 60 characters in length. Following this line is the NAMELIST 
data. The following is an abbreviated list of variable definitions. The 
NNEPEQ manual should be consulted for the complete list. 

NCOMP = The total number of components including controls that will be 
configured. 

TABLES = TRUE if maps are used, FALSE if not. (Default is TRUE) 

PUNT = Set PUNT = TRUE to use last good point as set of first guesses for 
next point. It is advisable to always have PUNT = TRUE. (Default 
is TRUE) 

LONG = Control for printing a history of the convergence process. It is 

advisable to have LONG = TRUE for new problems. (Default is TRUE) 

The first NAMELIST read process inputs these variables and then, if 
TABLES = TRUE, the code will read in the maps. Now the program is ready to 
read in the engine configuration data. 

Each type of component has a different set of input variables. However, 
the form of the input is the same for all components. For all components 

except controls the input data is read in the following form. 

KONFIG(l ,N)=M, JM1 , JM2, JP1 , JP2, 

SPEC(1,N)=V1,V2, .'. . , VI 5 

where N is the component number 
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M is 


JMl is 
JM2 is 
JP1 is 
JP2 is 


a number corresponding to the type of component where 

1 = inlet 

2 = duct or burner 

3 = gas generator 

4 = compressor 

5 = turbine 

6 = heat exchanger 

7 = splitter 

8 = mixer 

9 = nozzle 

10 = load 

11 = shaft 

12 = control 

the primary upstream airflow station number for flow components 
or the first component hooked onto a shaft 

the secondary upstream station number, or the second component 
hooked onto a shaft 

the primary downstream station number, or the third component 
hooked onto a shaft 

the secondary downstream station number, or the fourth component 
hooked onto a shaft 


The SPEC array is used to define the characteristics of the component 
identified in the KONFIG statement. The SPEC statement could specify as many 
as 15 different parameters but most components do not use all 15. Table 4.1 
is extracted from reference [ 3 ] and defines the values of the SPEC array for 
each component. 


Table 4.1 Definitions of the Conf igurat ion Variables 

Device Type 1: Inlet 

SPECCl ,N)=Inlet weight flow, lb/sec 
SPEC(2,N)=Free stream temperature, R 
SPEC(3,N)=Free stream pressure, psia 

SPEC(4 ,N)=Inlet drag table reference number, if blank computed 
SPEC(5,N)=Mach number at inlet 

SPEC(6 ,N)=Inlet recovery factor, constant or table reference number 
If =0, Mil spec is used 

SPEC(7 ,N)=Maximum permitted flow if table specified in SPEC(6,N) 
SPEC(8 ,N)=Scale factor on flow if table specified in SPEC(6,N) 

SPEC(9 ,N)=Altitude , feet, used only if SPEC(2,N), SPEC(3,N) =0 

SPEC( 10, N)= Fuel/ air ratio at inlet, usually =0 

SPEC( 11 ,N)=If nonzero, SPEC(9,N) is geopotential altitude 

SPEC( 1 2 , N)=Temperature change to be added to SPEC(2,N), usually zero 

SPEC(13,N)=Blank . 

SPEC( 14,N)=Corrected flow at inlet entrance (if specified, will 
override SPEC(l,N)) 

SPEC(l5,N)=Corrected flow at inlet exit (if specified, will override 
SPEC(1 ,N)) 

Device Type 2: Duct (or Diesel Engine) 

SPEC(1 ,N) -Pressure drop/inlet pressure or table reference number 
SPEC(2,N)=Optional , design duct Mach number 
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SPEC(3,N)=0 for duct, if nonzero then program assumes this device 

is a diesel engine and this is the mass of fuel injected per 
cylinder per cycle 

SPEC(4,N)=If SPEC(3,N)=0, then burner outlet temperature if positive 

simple duct if zero 
fuel/oxidizer ratio if negative 
If SPEC( 3 , N) . GT . 0 , then equal to diesel engine speed 
SPEC(5,N)=Burner efficiency or table reference number 
SPEC(6,N)=Fuel heating value or table reference number 
SPEC (7 ,N)=Cross sectional area of duct or. burner, sq inches 
SPEC(8 , N)=Rat io of inlet entrance bleed flow to total bleed available 
SPEC(9,N)=Exit bleed/total flow 
SPEC(10,N)=Fraction of air not heated 


Device Type 3: Gas Generator 

SPEC ( 1 , N)=Output temperature, R 
SPEC ( 2 ,N)=Fuel/oxidant weight ratio 
SPEC(3 ,N)=Generator pressure, psia 
SPEC(4,N)=Fuel flow rate, lb/sec 
SPEC(5,N)=Oxidant flow rate, lb/sec 

Not all of these SPEC's can be specified simultaneously. The options 

1. SPEC ( 1 ) and SPEC(4) 

2. SPEC(l) and SPEC(5) 

3. SPEC( 2) and SPEC(4) 

4. SPEC( 2) and SPEC(5) 

5. SPEC( 4 ) and SPEC(5) 

Device Type 4: Compressor 

SPEC(1,N)=R value used to read tables 
SPEC(2,N)=Compressor bleed flow/total flow 
SPEC(3,N)=Scale factor for corrected speed 
SPEC(4 ,N)=Table reference number for corrected speed 
SPEC(5,N)=Scale factor for corrected flow 
SPEC(6,N)=Table reference number for corrected flow 
SPEC (7 ,N)=Scale factor for efficiency 
SPEC(8 , N)=Table reference number for efficiency 
SPEC(9,N)=Scale factor on input pressure ratio, SPEC(l3,N) 
SPEC(10,N)=Third dimension argument value on map (variable geometr y) 
SPEC( 1 1 ,N)=Fractional bleed horsepower loss due to interstage bleed 
=0 means all bleed after full compression 
SPEC( 1 2 , N)=Des ired adiabatic efficiency at design point 
SPEC(l3,N)=Desired pressure ratio at design point R value and 
corrected speed 

SPEC( 14 ,N)=Corrected speed for design point on maps 


5: Turbine 

Pressure ratio at design point on maps 
Total bleed in/Totai bleed available 
Scale factor for corrected speed, usually 1 
: Table reference number for corrected flow 
: Scale factor for corrected flow 

: Table reference number for turbine adiabatic efficiency 
=Scale factor to get design point efficiency at design 
point on map 

SPEC(8 , N)=Scale factor calculated to get desired pressure ratio 


Device Type 
SPEC(1 ,N) 
SPEC( 2 , N) 
SPEC(3,N) 
SPEC(4,N) 
SPEC( 5 , N) 
SPEC( 6 , N) 
SPEC( 7 , N ) 


are 
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on map 

SPEC(9 ,N)=Turbine bleed flow at entrance/bleed flow 
SPEC(10,N)=Third dimension argument value on map 
SPEC( 1 1 ,N)=Desired efficiency at design point 
SPEC( 12,N)=Corrected speed at design point on map 
SPEC(l3,N)=Turbine horsepower split (usually =1) 

Device Type 6: Heat Exchanger 

SPEC(1 ,N)=Pressure Drop/lnlet pressure for main flow 
SPEC(2,N)=Pressure drop/lnlet pressure for secondary flow 
SPEC(3,N)=Temperature rise for main flow (main flow must be heated) 
SPEC(4,N)=Ef f ectiveness 

SPEC(5,N)=Scale factor on effectiveness 
Device Type 7: Splitter 

SPEC( 1 ,N)=Bypass ratio (Bypass flow/main flow) 

SPEC(2,N)=Pressure drop/inlet pressure for main stream 
SPEC(3,N)=Pressure drop/inlet pressure for bypass stream 

Device Type 8: Mixer 

SPEC( 1 ,N)=Inlet area of main flow calculated at design point 
SPEC(2,N)=Inlet area of secondary flow calculated at design point, 
if SPEC(7 , N)=2 , then throat area 

SPEC(3,N)=Total to static pressure ratio at main flow inlet if greater 
than 1, if less than 1 = Mach number at design point 
SPEC(4 , N)=veloc ity coefficient on mixed flow velocity, if SPEC(7,N)=2 
then momentum coefficient for ejector primary (l=ideal, less 
then l=less than ideal) 

SPEC(5,N)=If =1 total inlet area is held fixed as second area varies, 
if =0, independent areas 

SPEC(6,N)=This parameter can have a value of 0 or 1 . If =0 the definition 
of SPEC( 3 , N) is as described above. If =1 then SPEC(3,N) is the 
minimum total to static pressure ratio or the minimum Mach 
number of either stream 

SPEC(7 ,N)=I f =0 then normal mixer, =1 supersonic mixer, =2 ejector 

To use SPEC(7,N)=2, the user should consult the NNEPEQ manual 
for additional input requirements 

Device Type 9: Nozzle 

SPEC( 1 , N)=Flow area, sq inches, exit for converging, throat for C~D 
SPEC( 2 , N)=Flow coefficient or table reference number 
SPEC(3,N)=Blank 

SPEC(4 , N)=Nozzle exit static pressure or component no. (see SPEC(9,N)) 
SPEC(5,N)=Cv, velocity coefficient or table reference number 
SPEC(6,N)=0 if converging, =1 if C-D 

SPEC(7 ,N)=Area switch, =0 to fix area to input value, =1 to vary area 
to match flow required 
SPEC(8,N)=Blank 

SPEC(9,N)=If SPEC(4,N)=0, set SPEC(9,N) to component number of inlet 
Device Type 10: Load 

SPEC( 1 ,N)=load horsepower (negative) or table reference number 
SPEC ( 2 ,N)=propeller efficiency or =0 
SPEC(3,N)=thrust/SHP at SLS 
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Device Type 11: Shaft 

SPEC( 1 ,N)=Actual shaft RPM 

SPEC( 2 ,N)=Gear ratio, first component on shaft. Component rpm/shaft rpm 
SPEC( 3 , N)=Gear ratio, second component on shaft. Component rpm/shaft rpm 
SPEC(4,N)=Gear ratio, third component on shaft. Component rpm/shaft rpm 
SPEC( 5 , N)=Gear ratio, fourth component on shaft. Component rpm/shaft rpm 
SPEC(6,N)=Mechanical efficiency of first component 
SPEC (7 ,N)=Mechanical efficiency of second component 
SPEC (8 ,N)=Mechanical efficiency of third component 
SPEC (9 , N) =Me chan i. cal efficiency of fourth component 

Device Type 12: Control 

SPCNTL(1 ,N)=M for the SPEC(M,N) of device N which is to be varied^ 

SPCNTL( 2 , N)=N for the component number of the component being varied 
SPCNTL( 3 , N) = l 00 if a station property is to be varied 
=200 if a DATOUT is to be varied 
=400 if a performance property is to be varied 
SPCNTL(4,N)=Number of station property to be varied 
=1 weight flow 
=2 total pressure 
=3 total temperature 
=4 fuel to air ratio 
=5 corrected flow 
=6 Mach number 
=7 static pressure 
=8 interface corrected flow error 
or DATOUT(L) (see output section for values of L for 
each component) 

or performance property to be varied 
=1 total engine airflow 
=2 gross jet thrust 
=3 fuel flow 
=4 net jet thrust 
=5 TSFC 

=6 net thrust/airflow 

=7 total inlet drag 

=8 total brake shaft HP 

=9 net thrust with installation drags 

=10 net SFC 

=11 inlet drag (lip and spillage) 

=12 boattail drag 

SPCNTL( 5 , N)=f low station number if SPCNTL( 3 , N) = 100 
component number if SPCNTL( 3 , N) = 200 
=0 if SPCNTL( 3 , N)=400 
SPCNTL(6 ,N)=value to be achieved 

SPCNTL( 7 ,N)=tolerance as fraction of value, if =1 default value of 
0.001 will be used, if -0 control is turned off 
SPCNTL(8 ,N)=minimum allowable value, if zero ignored 
SPCNTL( 9 , N)=maximum allowable value, if zero ignored 


Figure 4.1 shows a sample turbomachinery input file. This is the file 
that corresponds to the V-6 Compound-Cycle engine discussed later in the first 
sample case. The only difference between this file and a conventional 
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TEST CASE: COMBINED CYCLE ENGINE 

&D ICEC=0 ,NC0DE=1 ,LONG=T,NMODES = l , SEND 
£D MQDE = 1 , 

KONFI G(1 ,1) =1,1, 0,2,0 ,SPEC(1 ,1)=7.S,4*0 ,1 ,5*0 ,0,0, 

KONFI Gd ,2) =4, 2, 0,3,0 ,SPEC(1, 2) =1 .8,0,1,1001,1,1002,1 ,1003,1 , 
2. ,0 ,. 35, 1 .05,1 , 

KONFI G(1 ,3) =7, 3,0 , 1 2 , 4 , SPECC1 , 3) =0 . 455 , 0 , 0 , 

KONFI G( 1 , 4) =4 ,4,0 , 5 , 0 , SPECC1 , 4) =1 . 3 , 0 , 1 , 1004 , . 5 , 1 005 , 1 , 1005 , 1 , 
0 0 . 79 9 • 9132 1 

KONFI G(1 ,5) =6, 12, 5, 13, 6, SPEC ( 1 , 5) = . 02 , .02,250. ,0.4,1, 

KONFI G(1 ,6) =2, 6,0 ,7,0 ,SPEC(1,S)=0 .1 ,0 ,0 .000145,5122. ,1,13300. , 
FARRAY ( 1 ,6) =2,1 

KONFI G(1 ,7) =5, 7, 0,3,0 , SPEC ( 1 , 7 ) =4 . 3 , 0 , 1 , 1 007 , 1 , 1 003 , 1 , 1 , 

0,1,. 345 ,5000,1, 

KONFI G(1 ,8) =2, 8,0 ,3,0 , SPEC(1 ,8)=0 ,0,0,0, .39,18300 . , 

FARRAY ( 1 ,8) =2,1 

KONFI G<1 ,3) =5, 3,0 ,10 ,0 ,SPEC(1 ,5) =1.8, 0,1, 10 03, 1,1010, 1,1, 
0,1,-348,5000,1, 

KONFI G(1 ,10) =3, 10, 0,11,0 ,SPEC(1 ,10)=0.,1,0.,0., .3850 ,1 . , 0 . , 

0. ,1. , 

KONFI Gd ,12) =3, 13,0 ,14,0 , SPEC (1 ,12) =0,1, 0,0. ,0.3850,1 . ,1 . , 

0. ,1. , 

KONFI Qd ,13) =11 ,2, 4, 7,0 , SPECC1 ,13) =3*1 . 

KONFI G(1 ,14) =11 ,3,15,0 ,0 ,SPEC(1 ,14) =3*1 . 

KONFI Gd ,15) =10 ,0 ,0,0 ,0 ,SPEC(1 ,15) =-200 ,0,0, 

KONFI G (1 , 1£) =1 2 , SPCNTLC 1 ,16) =1,5, 10 0,3, 10, 0,0 » 

KONFI G(1 , 17) =12 , SPCNTLC 1 ,17) =1,7, 100, 3, 3, 0,0, 

KONFI Gd , 1S)=12,SF'CNTL(1 ,18)=1 ,4,10 0 ,8,7,0 ,0 , 

KONF I G d , 20 ) =1 2 , SPCNTLC 1 , 20 ) =1 , 2 , 1 00 , 3 , 4 , 0 , 0 , 

KONFI G(1 , 21) =12, SPCNTLd ,21)=1 ,1 ,100,8,2,0 , 0 , 

KONFI G(1 ,22) =12, SPCNTLd , 22) =1 , 13 , 200 , 8 , 13 , 0 , 0 , 

KONFI G(1 , 23) =12, SPCNTLd , 23) =1 , 1 5 , 20 0 , 8 , 1 4 , 0 , 0 , 

KONFI G(1 , 24) =12, SPCNTLd , 24) =3 , 5 , 200 , 3 , 5 , 0 , 1 , 

KONFI G Cl ,26) =12, SPCNTLd , 26) =10 , 2 , 200 , 5 , 2 , 4 . 0 , 0 , 

SEND 

&D SPEC (9,1 6) =1 , SPEC (3, 17) =1 , SPECC 3 , 13) =1 , 

SPEC ( 9 , 20 ) =1 , SPECC 9 , 21 ) =1 , SPECC 9 , 22) =1 , SPECC 9 , 23) =1 , 

SPECC 3, 26) =1 , SEND 
&D SPECC 3,, 6) =0 .00012 SEND 
&D SPEC(3,6)=G. 00010 SEND 
&D SPEC(3,6)=0 .00003 SEND 
&D SPEC(3,6)=0 .00003 SEND 
&D SPEC(3,6)=0 .000075 SEND 
&D SPEC(3,6)=0 .00007 SEND 
&D SPECC3,6)=0 .000065 SEND 
&D SPECC3,6)=0 .00006 SEND 
&D SPEC(3,6)=0 .00005 SEND 
&D SPECC3,6)=0. 00004 SEND 
&D SPEC(3,6)=0 .00003 SEND 


Figure 4.1 Sample Turbomachinery Input File 
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turbomachinery input file for NNEPEQ is in the line defining the 
characteristics of component 6, repeated below. 

KONFIG( 1 , 6)=2 ,6 ,0, 7 ,0, SPEC(l ,6)=0. 1 ,0,0,000145 , 6122. , 1 , 18300. , 

This line appears to define the characteristics of a duct configured as a 
burner. However, the third item in the SPEC list for a duct is listed as 
unused by the NNEPEQ manual [3, p.6l. In this case, this variable has been 
assigned to be an indicator of the presence of a diesel engine in the engine 
configuration. If the value of this variable is zero, then the component is 
assumed to be a conventional burner. If the value is non-zero then the 
component is assumed to be a diesel engine and the numerical value of this 
variable corresponds to the specified value of fuel injected into each 
cylinder per cycle (lb/cycle). Further, if the value of the third variable is 
non-zero, then the value of the fourth variable, ordinarily assigned to be the 
desired burner outlet temperature, becomes the specified value of diesel 
engine speed. When the value of the third variable is non-zero, the 
other variables in the SPEC array are not used. 

When developing a data set for a new engine configuration, it is 
suggested that the third variable be set equal to zero and estimated values of 
the diesel engine pressure drop and exhaust temperature be assigned. This 
will allow the turbomachinery portion of the program to be tested without the 
added complication of the diesel program. 

Figure 4.2 shows the form in which data must be entered for the diesel 
data set. The data is read format-free but must be in the order shown. The 
variables in the input data set are defined as shown in Table 4.2. 


Table 4.2 Definitions of Variables in the Diesel Input File 


TITLE 

NCYCLE 

NCYL 

BORE 

STROKE 

CONROD 

CR 

MODE 

Cl 

C2 

C3 

C4 

C5 

ANNND 

THEAD 

TPISTN 

TSLEEV 

MEXH 

EVO 

EVC 

MINT 


Title of run, up to 80 characters 
2 if two stroke, = 4 if four stroke engine 
number of cylinders 
Cylinder bore, (inches) 

Stroke, (inches) 

Connecting rod length, (inches) 

Geometric compression ratio. 

Piston-cylinder mode. ( = 1 for conventional piston-cylinder 
configuration, =2 for opposed piston configuration) 

Shape parameter for diffusion burning. 

Crankangle for start of combustion. 

Combustion duration parameter. 

Premixed burning fraction. 

Shape parameter for premixed burning. 

Multiplier for Annand heat transfer correlation. 

Cylinder head temperature, (deg R) 

Piston temperature, (deg R) 

Cylinder wall temperature, (deg R) 

Exhaust valve or port indicator. (=0 for exhaust ports 
and =1 for exhaust valves) 

Exhaust valve or port opening crankangle. (degrees ATDC) 
Exhaust valve or port closing crankangle. (degrees ATDC) 
Intake valve or port indicator. (=0 for intake ports 
and =1 for intake valves) 
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TITLE 
N CYCLE 

NCYL 




BORE 

STROKE 

CON ROD 

CR ’ 

MODE 

Cl 

ANNND 

C2 

THEAD 

C3 

TP I STM 

C4 

TSLEEU 

C5 

MEXH 
NT EXH 
ALPHEXi 
ALPH EX' 

ECO 
CD EXH 

;i) FEXH(l) 
:;2) FEXH( 2) 

ECO 

M I NT 

ACO 


ALPH EX < NT EX H ) F EXH < NT EX H ) 
WIDTH! CD I NT 


Figure 4.2 Format of Diesel Input File 



AVO 

AVC 

NTEXH 

CDEXH 

ALPHEX(I) 

FEXH(l) 

WIDTHE 

NTINT 

CDINT 

ALPHIN(I) 

FINT(I) 

WIDTHI 


Intake valve or port opening crankangle. (degrees ATDC) 

Intake valve or port closing crankangle. (degrees ATDC) 

Number of exhaust valve area versus crankangle data pairs 
to be read in. (for MEXH=1 only) 

Exhaust valve or port discharge coefficient 

Array of crankangles for exhaust valve flow areas, (deg ATDC) 
Array of exhaust valve flow areas, (square inches) 

Fraction of the cylinder circumference devoted to exhaust port. 
Number of intake valve area versus crankangle data pairs 
to be read in. (for MINT=1 only) 

Intake valve or port discharge coefficient 

Array of crankangles for intake valve flow areas, (deg ATDC ) 
Array of intake valve flow areas, (square inches) 

Fraction of the cylinder circumference devoted to intake port. 


When MEXH=1, exhaust valves are being used and the used must specify 
NTEXH, CDEXH and the arrays of crankangle versus flow area, ALPHEX(l) and 
FEXH(l). When MEXH=0, exhaust ports are used and the user need only specify 
WIDTHE and CDEXH. ALPHEX(l) and FEXH(l) are calculated internally. A 
similar procedure is used when MINT=0. The case illustrated in Table 1 
corresponds to an engine that has exhaust valves and intake ports. 
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5 Output from the Combined-Cycle Engine Simulation Program 

The output from the combined-cycle engine simulation program consists of 
two files. The first file characterizes the turbomachinery portion of the 
engine and is called TURBO. OUT and the second file characterizes the diesel 
engine and it is called DSL. OUT. The first file is identical to the output 
file from NNEPEQ [2,3]. The input and output station properties of the 
diesel engine are reported in this file as if the engine was a burner. This 
file contains a listing of the properties such as mass flow rate, temperature, 
pressure, fuel/air ratio, etc. at each flow station for the engine. Then the 
output variables that are individually defined for each component are 
presented in a table. The table identifies these variables by their position 
in the DATOUT array calculated by CCE. The definitions of the elements of 
this array are provided in Table 5.1. Finally, the file includes a summary of 
the engine as a whole including air flow rate, fuel consumption, thrust, etc. 


Table 5.1 Definition of the elements of the DATOUT array 

Device Type 1: Inlet 
DATOUT( l)=inlet drag 
DATOUT(2)=velocity, ft/sec 
DATOUT(3)=velocity, knots 
DATOUT(4)=ram temperature ratio 
DATOUT(5)=ram pressure ratio 
DATOUT(6)=Mach number 

DATOUT(7)=Inlet recovery - exit total pressure/ram pressure 
DATOUT(8)=exit temperature/ 518.67 
DATOUT(9)=altitude, feet 

Device Type 2: Duct or Diesel engine (Diesel engine output is provided 
in a separate data file) 

DATOUT( 1 )=Pressure drop/inlet pressure from momentum 
DATOUT( 2)=Pres sur e drop/inlet pressure from SPEC(l,N) 
DATOUT(3)=Pressure ratio at duct inlet used to compute inlet Mach 
number (Total/Static) 

DATOUT ( 4 ) = Fue 1 flow/duct inlet weight flow 
DATOUTC 5)=Cross sectional area - sq inches 
DATOUT ( 6 )=Fuel flow, lb/hr 

DATOUT( 7 )=In let Mach number (if SPEC(2,N) or SPEC(7,N) were specified 
at the design point) 

DATOUT(8)=Burner efficiency 

DATOUT(9)=Burner outlet temperature (before bypass added) 

Device Type 3: Gas Generator 

DATOUT( 1 )=Generator temperature, deg R 
DATOUT(2)=As signed generator fuel/oxidizer ratio 
DATOUT( 3)=Generator pressure, psia 
DATOUT(4)=Generator fuel weight flow, lb/sec 
DATOUT( 5)=Generator oxidant weight flow, lb/sec 
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DATOUT(6)=Total generator weight flow, lb/hr 
DATOUT (7)=Calculated generator fuel/oxidizer ratio 

Device Type A: Compressor 

DATOUT( 1 )=Hor sepower required (negative) 

DATOUTC 2) =Phys i cal RPM 

DATOUT(3)=Third dimensional argument on compressor maps 

DATOUT(4)=R value used on maps 

DATOUT(5)=Surge margin in percent 

DATOUT(6)=Corrected speed used to read maps 

DATOUT(7)=Scale factor on corrected flow 

DATOUT (8)=Compressor efficiency 

DATOUT(9)=Compressor pressure ratio 

Device Type 5: Turbine 

DATOUTC 1 )=hor sepower produced by turbine 
DATOUT( 2)=phys ica 1 rpm 

DATOUT( 3)=Third dimension argument value on turbine maps 

DATOUT(4)=pressure ratio used in table lookup 

DATOUTC 5)=scale factor on corrected speed 

DATOUT(6')=Corrected speed used in table lookup 

DATOUTC 7 )=scale factor on corrected flow 

DATOUTC 8) “turbine efficiency 

DATOUTC 9) “turbine overall pressure ratio 

Device Type 6: Heat Exchanger 

DATOUTC l)=Pressure drop/inlet pressure for main flow 
DATOUT(2)=Pres sure drop/inlet pressure for secondary flow 
DATOUTC 3 )=B lank 
DATOUTC 4 )=Ef fectiveness 

DATOUTC 5 )=S cal e factor on effectiveness 
DATOUTC 6) “Calculated temperature difference 
DATOUTC 7) “Temperature dif ference/(T hot - T cold) 

DATOUT(8)=Temperature rise difference (guess value/calculated - 1) 

Device Type 7: Splitter 
DATOUTC l)=Bypass ratio 

DATOUT(2)=Pressure drop/inlet pressure in primary flow stream 
DATOUTC 3 )=Pres sure drop/inlet pressure in secondary flow stream 

Device Type 8: Mixer 

DATOUTC 1 )=Ma in flow area - sq inches 
DATOUT(2)=Secondary flow area - sq inches 

DATOUTC 3) “Total to static pressure ratio at main flow inlet 
DATOUT(4)=Total to static pressure ratio at secondary flow inlet 
DATOUTC 5 )=Veloc ity at main flow inlet 
DATOUTC 6)“Ve loc ity at secondary flow inlet 
DATOUTC 7) =Ex it mixed flow velocity 

DATOUT(8)=Static pressure difference between streams 
DATOUTC 9 )=Total mixed to average static pressure ratio 

Device Type 9: Nozzle 

DATOUTC 1 )=Gross jet thrust, lb 

DATOUTC 2) “Actual jet velocity, ft/sec 

DATOUTC 3) “Total to static pressure ratio at throat 
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DATOUT(4)=Nozzle exit area, sq inches 
DATOUT(5)=Nozzle throat area, sq inches 
DAT0UT(6)=Cd , flow coefficient 
DATOUT(7)=Cv, velocity coefficient 
DATOUT (8)=Critical pressure ratio at throat 
DATOUT(9)=Overal 1 pressure ratio, inlet total 


to exit static 


Device Type 10: Load 

DATOUT( 1 )=Load horsepower (negative) 
DATOUT( 2)=Actual shaft rpm 
DATOUT( 3)=Prope ller thrust 


Device Type 11: Shaft 
DATOUT ( 1 )=Net shaft ho 
DATOUT (2)=Actual shaft 
DATOUT( 3)=Ac tual shaft 
DATOUT(4)=Actual shaft 
DATOUT( 5)=Actual shaft 
DATOUT ( 6 ) = Ac t ua 1 shaft 
DATOUT(7)=Blank 
DATOUT(8)=Net shaft ho 

Device Type 12: Control 
There is no DATOUT arr 


sepower (required-delivered) 
rpm 

rpm of first component 
rpm of second component 
rpm of third component 
rpm of fourth component 

sepower/total horsepower 
y for controls 


The second output file shows the calculated results for the diesel 
engine. It consists of two parts and is essentially identical to the output 
file from DSL2 [ 1 ] . The first part is a summary of the calculated quantities 
that characterize the diesel engine's performance. The second part is a 
listing of the cylinder pressure and temperature and other quantities at 
10 degree intervals during the engine cycle. Both parts for the diesel 
engine output listing are self-explanatory although the symbols are defined in 
Table 5.2 below. The only difference between the diesel output file from CCE 
and the output file from DSL2 is that DSL2 gives information about only a 
single cylinder of a multi-cylinder engine while CCE gives information about 
the entire engine. 


Table 5.2 Symbols Used in Diesel Output File 

IMEP = Indicated Mean Effective Pressure 
BMEP = Brake Mean Effective Pressure 
FMEP = Friction Mean Effective Pressure 
I HP = Indicated Horsepower 
BHP = Brake Horsepower 
FHP = Friction Horsepower 

ISFC = Indicated Specific Fuel Consumption 
BSFC = Brake Specific Fuel Consumption 

EQUIVALENCE RATIO DURING COMPRESSION = This is the equivalence ratio 

in the cylinder at the start of the compression process. It is 
an indicator of the amount of residual gas left from the previous 


18 



cycle. 

EQUIVALENCE RATIO AT EVO = This is the equivalence ratio in the cylinder 
when the exhaust valve (or port) opens. It is an indicator of how 
much of the oxygen that entered the cylinder is still unused. 

EQUIVLANCE RATIO BASED ON FUEL AND AIR FLOW RATES = This is the equivalence 
ratio calculated from the fuel and air flow rates supplied to the 
engine. In two _ stroke engines it will generally be less than that 
in the cylinder at EVO due to the diluting effect of the extra air 
supplied to scavenge the exhaust products from the cylinder. 


CA 

PCYL = 
TCYL = 
PHICYL = 
MCYL = 
MDOTIN - 
MDOTOUT= 
BALANCE= 


Crankangle 
Cylinder pressure 
Cylinder temperature 
Equivalence ratio in the cylinder 
Total mass in the cylinder 

Mass flow rate into cylinder through intake valve (or port) 

Mass flow rate out of cylinder through exhaust valve (or port) 
Monitor on the overall cylinder energy balance. (Generally less 
than 0.005) 


19 



6 Sample Cases 


6.1 Case Number One: V-6 Diesel Engine With Aftercooling 

This case is derived from the combined-cycle engine proposed by Castor 
of the Garrett Turbine Engine Company [7]. This engine assumes a 6-cylinder 
diesel engine running at 6122 rpm coupled to turbomachinery that provides 
supercharging air at 150 psia to the diesel as well as power recovery from the 
exhaust. The overall engine power is designed to be 1000 hp. The engine 
specifications are shown in the table below. 


Table 6.1 V-6 Engine Specifications 


Two-stroke, V-6 Diesel Engine, uniflow 
Bore 3.10 inches 

Stroke 2.94 inches 

Compression ratio 9.17 (geometric) 

At design point (14.7 psia, 519 deg R 
Inlet manifold pressure 

Inlet manifold temperature 

Compressor efficiency 
Turbine efficiencies 
Exhaust opens 
Exhaust closes 
Intake opens 
Intake closes 


scavenged 


intake conditions 
150 psia 
896 deg R 
0,79 

0.845, 0.848 
90 deg ATDC 
239 deg ATDC 
126 deg ATDC 
234 deg ATDC 


) 


The block diagram shown in Figure 6.1 shows the combined-cycle engine 
configuration used to model the engine. The numbers enclosed by circles 
identify "flow stations" in the engine and the numbers and names in the 
rectangles identify components. The design includes a low pressure fan to 
supply air through a splitter to both the high pressure compressor and an 
air-to-air aftercooler. The diesel engine receives air from the aftercooler 
and discharges to a turbine that drives the fan and compressor. Exhaust 
products pass from this turbine through a duct and into a power turbine. The 
duct was included so that reheating of the exhaust products before entering 
the power turbine could be investigated. However, this feature is not used 
here. 

The engine simulation program was adjusted to match the performance of 
the engine proposed by Castor [7] at the design point. Then, with the diesel 
speed held constant, the fuel flow rate was reduced to provide information 
about part load performance. Figure 6.2 shows the turbomachinery input data 
file used by CCE. The corresponding diesel input file is shown in Figure 6.3 

Table 6.2 shows a comparison between the performance characteristics 
listed by Castor [7] and those provided by CCE. 


20 



t c> 

















TEST CASE: COMBINED CYCLE ENGINE 

£D I CEC=0 ,NCODE=l , LGNG=T ,NMODES=l , &END 
6D MODE=l , 

KONFIG(l,l)=l ,1 ,0 ,2,0 ,SPEC(1 ,1)«7.8,4*0 ,1 ,5*0 ,0 , 0 , 

KONFI G( 1 , 2) =4 , 2 , 0 , 3 ,0 , SPEC Cl , 2) =1 .8,0 ,1 ,1001 ,1 ,1002 ,1 ,1003,1 , 
2. ,0 , .85,1 .05,1 , 

KONFI G(1 ,3)=7 >3,0 ,12 , 4,SPEC(1 , 3) =0 .455,0 ,0 , 

KONFI G(1 ,4) =4, 4,0 , 5 , 0 , SPEC( 1 , 4) =1 . 3 , 0 , 1 , 1004 , .5,1005,1,1006,1, 
0,0, .73,9.9132,1, 

KONFI G( 1 , 5) =6 ,12,5,13,6, SPEC( 1 , 5) = . 02 , .02,250. ,0.4,1, 

KONFI G(1 ,6) =2, 6,0 ,7,0 , SPEC( 1 , 6) =0 . 1 , 0 , 0 . 000145 , 61 22 . ,1 ,18300 . , 
FARRAY(1 ,6) =2,1 

KONFI G(1 ,7) =5, 7,0 ,B, 0 , SPECC 1 , 7) =4 . 3 , 0 , 1 , 1007 , 1 ,1008,1,1, 

0,1,. 845,5000,1 , 

KONFI G(1 ,S)=2,8,0 ,9,0 , SPECC 1 ,S)=0 ,0 ,0 ,0 , .33,18300 . , 

FARRAY ( 1 ,8) =2,1 

KONFI G(1 ,3) =5 ,9, 0 ,10,0, SPECC 1 , 3) =1 . 8 , 0 , 1 , 1003 , 1 ,1010 ,1,1, 

0,1, .848,5000,1, 

KONFI G(1 ,10) =3, 10 ,0,11,0, SPEC Cl , 10 ) =0 . , 1 , 0 . , 0 . , .3850 , 1 . , 0 . , 

0 . , 1 . , 

KONFI G(1 ,12) =3, 13, 0,14,0 , SPECC 1 , 12) =0 , 1 , 0 , 0 . ,0.3850 ,1 . , 1 . , 

0 . , 1 . , 

KONFI 6( 1,13) =11, 2, 4, 7,0, SPEC ( 1 , 13) =5*1 . 

KONFI G(l, 14) =11, 9, 15, 0,0, SF'EC( 1,14) =3*1 . 

KONFI G(1 ,15) =10 ,0 ,0 ,0 , 0 ,SPEC(1 ,15) =-200 ,0,0, 

KONFI G( 1 , 16) =12, SPCNTL ci , 1 6) =1 , 3 , 1 00 , 8 , 1 0 , 0 , 0 , 

KONFIGC1 ,17)=12,SPCNTL(1 ,17)=1 ,7,100 ,8,3,6 ,0 , 

KONFI G(1 , 1S)=12,SF'CNTL( 1 ,18)=1 ,4,100 ,8,7,0 ,0 , 

KONFI 6(1 , 20 ) =12 , SPCNTLC 1 , 20 ) =1 , 2 , 100 , 8 , 4 , 0 , 0 , 

KONFI G(1 , 21 ) =1 2 , SPCNTL C 1 , 21 ) =1 , 1 , 100 , 8 , 2 , 0 , 0 , 

KONFI G(1 , 22) =12, SPCNTL (1 , 22) =1 , 1 3 , 200 , 8 , 1 3 , 0 , 0 , 

KONFI GC1 , 23) =12 , SPCNTL ( 1 , 23) =1 ,15,200,8,14,0,0, 

KONFI GC1 ,24) =12, SPCNTL (1 ,24) =3, 5, 200 ,8,5,0 ,1 , 

KONFI G(1 , 26) =12, SPCNTL (1 ,26) =10 ,2,200 ,5, 2, 4.0 ,0 , 

SEND 

£D SPEC ( 9 , 1 6 ) =1 , SPECC 3 , 1 7 ) =1 , SPECC 9 , 18) =1 , 

SPEC ( 9 , 20 ) =1 , SPEC ( 9 , 21 ) =1 , SPECC 9 , 22) =1 , SPECC 9 , 23) =1 , 

SPECC 9 , 26) =1 , SEND 
&D SF'EC(3,6)=0 .00012 SEND 
&D SPEC(3,6)=0 .00010 &END 
&D SPEC(3,6)=0 .00003 SEND 
&D SPEC(3,6)=0 .00008 &END 
£D SPEC(3,6)=0 .000075 SEND 
£D SPEC(3,6)=0 .00007 Si END 
&D SPEC(3,6)=0 .000065 S.END 
&D SPEC(3,6)=0 .00006 S.END 
&D SPEC(3,6)=0 .00005 SEND 
£D SPEC( 3 , 6) =0 .00004 6 END 
&D SPEC(3,6)=0 .00003 SEND 


Figure 6.2 Turbomachinery Input Data File for Sample Case 1 
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TW 0 - ST P. OKE , U - 6 EN GINE N I T H UN I F L ON SCAM EN G I N G 
2 6 

3.100 2.340 7. 168 3.1712 1 

1 . 3 345 . 75 . 0.0 3.5 

0.55 1460. 1460. 1460. 

1 30. 233. 0 126. 234. 

13 0.8864 


30 . 00 

0 . 0000 

100 . 37 

0 .1343 

115.00 

0 . 5242 

121 . 61 

0 . 8514 

131.00 

1 .2346 

142.25 

1 . 5233 

1 62 . 83 

1 . 6324 

180 . 35 

1 .6000 

134.00 

1.3231 

201 . 55 

1 . 0338 

203.00 

0 . 6638 

222 . 23 

0 . 2617 

235.00 

0.0000 

0 .5542 

0 . 8 


Figure 6.3 Diesel Input Data File for Sample Case 1 
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Table 6.2 Combined-Cycle Engine Performance Comparison 



Castor [7] 

CCE 

Brake mean effective pressure (psi) 

393 

391.7 

Brake power (diesel only) (HP) 

810 

806.5 

Brake power (power turbine) (HP) 

241 

200.0 

Brake power (total) (HP) 

1051 

1006 . 5 

Brake specific fuel consumption (lb/hp-hr) 

0.328 

0.317 

Maximum cylinder pressure (psia) 

3362 

3140. 

Air flow rate (lb/sec) 

2.44 

2.44 

Turbine inlet temperature (deg R) 

2214 

2208. 

Equivalence ratio 

0.68 

0.674 


The diesel output file for the design point case is shown in Figure 6.4 and 
the turbomachinery output file is shown in Figure 6.5. The overall BSFC shown 
in Table 6.2 is calculated from the power produced by both and diesel and 
power turbine and the fuel consumed by the diesel. It is not calculated 
within the program CCE. 

The turbomachinery input file shown in Figure 6.2 specifies that the 
program will evaluate the design point case as well as 11 off-design cases. 
Although performance may only be desired at a few off-design points, in 
practice it is usually necessary to run the program at several intermediate 
points. The program uses an iterative procedure to solve the turbomachinery 
matching problem. The previous case conditions are used as initial estimates 
for the next case. If the change in conditions from one case to the next is 
too large, the program will be unable to converge. When this is a problem, 
the solution is usually as simple as inserting several intermediate points 
between the desired cases. 

In this case, the part load, rated speed performance was investigated. 
Table 6.3 shows the effect of reducing the amount of fuel injected on the 
power split between the diesel engine and the power turbine as well as the 
effect on the overall brake specific fuel consumption. 


Table 6.3 Effect of fuel flow rate on Combined-Cycle Engine Performance 


Fuel injected 
per cylinder per cycle 

BMEP 

Brake Power 
(diesel ) 

Brake Power 
(turbine) 

BSFC 

( lb/hp-hr) 

(lb) 

(psi) 

(HP) 

(HP) 

0.000145 

391.7 

806.5 

200.0 

0.317 

0.00012 

325.7 

670.6 

156.6 

0.320 

0.00010 

271.1 

558.2 

125.0 

0.322 

0.00008 

214.2 

441.0 

89.6 

0.332 

0.00006 

156.2 

321.5 

55.7 

0.351 

0.00004 

97.1 

200.0 

26.4 

0.389 
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TWO-STROKE, V-6 ENGINE WITH 

BORE 

STROKE 

CONNECTING ROD 
GEOMETRIC COMPRESSION RATIO 
EFFECTIVE COMPRESSION RATIO 
NUMBER OF CYLINDERS 
SWEPT VOLUME (TOTAL) 

ENGINE SPEED 
MEAN PISTON SPEED 
SUPPLY AIR PRESSURE 
SUPPLY AIR TEMPERATURE 
EXHAUST PRESSURE 


UN I FLOW SCAVENGING 
3.100 INCHES 


2.940 

7.168 

9.17 

7.50 


INCHES 

INCHES 


133.141 CUBIC INCHES 
6122. RPM 
3000. FT/MIN 
150.00 PSIA 
885.86 DEG R 
133.74 PSIA 


PORT TIMING: INTAKE OPEN 126. 

INTAKE CLOSE 234. 
EXHAUST OPEN 90 . 
EXHAUST CLOSE 239. 


FUEL FLOW RATE 

0.0888 LBM/SEC 

FUEL INJECTED/CYCLE 

0 . 000145 

LBM/ CYCLE 

AIR FLOW RATE 

2.4384 LBM/SEC 

AIR INDUCTED/ CYCLE 

0 . 003983 

LBM/ CYCLE 

IMEF' 

419.7 PSI 

I HP 

864.2 HP 


BMEP 

391 .7 PSI 

BHP 

806.5 HP 


FMEP 

28.0 PSI 

FHP 

57.7 HP 



ISFC 

BSFC 


0.3698 LBM/HP-HR 
0.3962 LBM/HP-HR 


INDICATED THERMAL EFFICIENCY 0.3779 

BRAKE THERMAL EFFICIENCY 0.3527 


MAXIMUM PRESSURE 
MAXIMUM TEMPERATURE 
MAXIMUM RATE OF PRESS RISE 


3140.0 PSIA AT 371.3 DEG 

4020.8 DEG R AT 390.0 DEG 

60.50 PSI/DEG AT 355.3 DEG 


PURITY 0.8371 
SCAVENGING EFFICIENCY 0.7759 
TRAPPING EFFICIENCY 0.7763 
CHARGING EFFICIENCY 0.5268 
DELIVERY RATIO 0.6787 
RESIDUAL FRACTION 0.2366 


EQUIVALENCE RATIO DURING COMPRESSION 0.1539 
EQUIVALENCE RATIO AT EVO 0.6738 
EQUIVALENCE RATIO BASED ON FUEL AND 

AIR FLOW RATES (EXHAUST) 0.5234 


% FUEL ENERGY TO HEAT LOSS 6.67 % 
% FUEL ENERGY TO WORK 37.79 % 
% FUEL ENERGY TO EXHAUST 55.54 % 

EXHAUST TEMPERATURE 2208.0 DEG R 


Figure 6.4 Diesel Output File for Sample Case 1 
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* 


CA 

PCYL 

90.00 

545.3 

100 .00 

453.6 

110.00 

377.1 

120.00 

292.3 

126.00 

241 .0 

130 . 00 

205.7 

140.00 

142.8 

150.00 

137.0 

160.00 

133.7 

170.00 

140.3 

175.17 

141 .0 

180.00 

142.9 

180.00 

142.9 

190.00 

144.9 

200.00 

146.8 

210 . 00 

149.3 

220.00 

150.4 

230 . 00 

153.9 

240 .00 

1 6 6 . 5 

250.00 

136.2 

260.00 

213.1 

270.00 

250 .2 

280.00 

302.2 

290.00 

376.5 

300.00 

434.3 

310 . 00 

645.6 

320.00 

386.0 

330.00 

1236.4 

340.00 

1693.9 

350 .00 

2193.3 

360 . 0 0 

2735.0 

370.00 

3133.9 

380.00 

2935.3 

390.00 

2337 . 7 

400.00 

1810.2 

410 . 00 

1353.3 

420.00 

1030.0 

430 . 00 

307.7 

440 .00 

653.3 


TCYL 

PHICYL 

3047. S 

0 .6738 

2930 .4 

0 .6733 

2303.5 

0.6733 

2647.3 

0.6733 

2534.3 

0.6733 

2445.0 

0 .6738 

2206.3 

0 . 6483 

1879.3 

0 .4754 

1641 .3 

0.3476 

1435.7 

0 .2632 

1427.9 

0.2395 

1333.3 

0 .2136 

1333.3 

0.2136 

1324.4 

0.1869 

1232.3 

0 .1662 

1263.9 

0.1553 

1263.9 

0.1539 

1272.1 

0.1539 

1299.2 

0.1539 

1333.1 

0 .1539 

1386.2 

0.1539 

1445.0 

0 .1539 

1516.5 

0.1539 

1603.1 

0.1539 

1707.4 

0 .1539 

1332.0 

0 .1539 

1978.0 

0 .1539 

2141 .6 

0.1539 

2307.1 

0.1539 

2476.1 

0 .1610 

2904.7 

0.2355 

34S2.7 

0 .3746 

3898.6 

0.5170 

4020 .3 

0.6125 

3918 . 3 

0 .6565 

3719.3 

0.6705 

3512.9 

0.6736 

3330.3 

0 .6741 

3176.9 

0 .6741 


MCYL 
0.004199 
0.004137 
0.003929 
0,003507 
0.003145 
0.001537 
0.001459 
0.001967 
0.002667 
0.003438 
0.003842 
0 .004010 
0.004011 
0.004239 
0.004360 
0.00437 7 
0.004239 
0.004088 
0.004051 
0.004051 
0 .004051 
0.004051 
0.004051 
0 .004051 
0.004051 
0.004051 
0.004051 
0 .004051 
0.004051 
0.004053 
0.004073 
0.004112 
0 .004152 
0,004173 
0.004191 
0.004195 
0.004195 
0.004196 
0.004196 


MOOT IN 
0 .000000 
0.000000 
0.000000 
0.000000 
0.000000 
-0.000016 
0.000031 
0 .000063 
0.000075 
0.000073 
0.000078 
0.000070 
0.000070 
0.000058 
0.000041 
0.000016 
-0.000006 
-0.000006 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0,000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0 .000000 
0 .000000 
0.000000 


MOOT OUT 
0.000000 
0 .000013 
0.000030 
0 .000055 
0.000064 
0.000063 
0 .000029 
0 .000019 
0.000024 
0.000028 
0 .000029 
0 .000040 
0 .000040 
0 .000040 
0.000034 
0 .000021 
0.000011 
0 .000005 
0.000000 
0 .000000 
0.000000 
0.000000 
0 .000000 
0.000000 
0 . 000000 
0 .000000 
0.000000 
0.000000 
0.000000 
0 .000000 
0.000000 
0 .000000 
0.00000 0 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 


BALANCE 
0 .00000 
0.00000 
0 .00000 
0.00000 
0 .00000 
0.00000 
0 .00000 
0.00000 
-0.00003 
- 0 .00001 
- 0.00002 
0.00006 
-0 .00003 
-0 . 00003 
-0 .00004 
-0.00005 
-0 .00004 
-0.00004 
-0 .00005 
-0 .00005 
-0 .00005 
-0.00005 
-0.00005 
-0 .00005 
-0.00005 
-0.00005 
-0 .00005 
-0 .00005 
-0.00005 
-0 .00005 
-0.00004 
-0 , 00004 
-0 .00005 
-0.00004 
-0.00004 
-0 .00004 
-0 .00004 
-0 .00004 
-0 .00005 


Figure 6.4 - - continued 
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CASE IDENTIFICATION TEST CASEs COMBINED CYCLE ENGINE 


rHOrHO^TrHOrHrHrHOOO 

OOOOOOOOOOOOO 

AAAAAAAA AAAA A 

oo(MooCDo«-iNoooo 

oolTioCOOoo^fftooo 

inoHOooo^^Nooo 

OO©O(M©oC0CMOOOO 

HOIJlOWWOHriHOOO 


OOOOOOOOOOOOOO 


o o 


LJ 

CJ 

2 

at 

LJ 


at 

O o 
at © 
at CD + 

3 3© 
H O 

S <L o 
I- O 
JCftO 
Li. O 


O O O 

o o o 

aaaa 

o o o - 
o o o 
o o o 
o o o 
o o o 


o o 
o o 
' + 
© 
© o 
o o 

O O 

o o 
o o 


££££££££ 


o o 
o o 
* + 
© 
o o 
o o 
o o 
o o 
o o 


o o 
o o 
+ 
© 
o o 
o o 
o O 
o o 
o o 


o o o o 
o o o o 
o o o o 
o o o o 
o o o o 


00 + 
1 - © 
3 o 
o o 
I- O 


£ £ 
o o 
o o 
o o 
m o 
CD O 


o in rt o o 
o «H O O O 

AAAAA 

o <\J o o o 
o o o o o 

O H O If) O 
© «H O <7 O 

N H H CO O 


O H H O 
o o o © 

AAA A 

O W H O 

o m d o 
CO ID © O 
<7 ffl 00 O 
00 H H O 


AA 

o o 
o o 
o o 
o o 
© © 


at c- 


OOOOOOOOOOOOOO 
OOOOOOOOONNOfJN 

oooooooooooooo 

=,£: AAAAA A AAAAAAAA 

*-«CD|-OOOOOOOOO‘Xi's0OCDL0 
HCD<looooooooo©©o©© 
^UhOOOOOOOOO^lCO'i'f) 
l-a:tf)ooociooooo9vo^^ 
to CL ooooooooorHrHCH^ 

oooooooooooooo 


oooooooooooooo 


N + 

O o 

t- O 

22 


AA 

ffl O 
ID O 
<7 O 
OJ o 
N O 


l 

rH O O rH 

o o o o 

A A A A 

t-i O O ffl 
ID CO O © 
CD H O o 
c\j o © 
WHOM 


o «H 
o o 

AA 

o © 

O © 
O CM 
o CM 
o CO 


o o 
O O 

AA 

O o 
o o 
ID ID 
CD CD 
© © 


o o o 
o o O 

AAA 

o o o 
o o o 
o o o 
o o o 
o o o 


<r 


5 

CL 

5 

o 

> 

h~ 

at 

LJ 

0. 

o 

Ct 


<r 


at 
U 
X CD 

° £ 

2i 


Ci 

LJ 

at 


L0 + • 
3 © i 


, A A A A AAA 


© 6 © © © o © 

t-ooooooooo»Hffl 
<Xoooooooooo© 

1- O 
CO O 


ooooooooCDfO 
ooooooooo© 


o o o 
o o o 

AAA 

O CD CTt 
o © CO 
o CM © 
O O © 
OWri 


oooooooooooooo 


O rt O H (O 00 V O ^ rH rH © O O 
OOOOOOOOOOOOOO 


2AAAAAAAAAAAAAA 

DOOOONNOOOOOOOO 


o 

O 

o 

O 

O 

O 

o 

o 

o 

o 

O 

o 


o 

O 

o 

o 



o 

o 

O 




f— ( 

rH 


H 


rH 

rH 

O 

CM 

o 

o 

O 

o 

O 

o 

O 

o 

o 

O 

o 

o 


O 

o 

o 

o 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


ID + 

+ 

+ 

+ 

o 

£ 

o 

£ 

© 

£ 

O 

£ 

© 

£ 

© 

© 


t O 

© 

o 

© 

^7 

*rt 

ID 

H 

O 

<7 

^7 


l£i 

<7 

<7 


<r 

3 o 

CM 

o 

CJ 

© 

CM 

CD 

CD 

CJ 

CJ 

CJ 

CD 

CD 

N 

«rt 

H 

H 

o o 

CD 

o 

CM 

CO 

^7 

in 

CM 

CD 

iO 

iO 

© 

© 

<7 

O 

o 


H o 

li) 

o 

CM 

<7 

CD 

«7 

rH 

N 

© 

© 

CD 

CD 

rH 

© 

© 

© 

2 ° 

ID 

o 

o 

IS 

CM 

CD 

CD 

ID 

«rt 

H 

n 

CD 

ID 

ID 

ID 

1- 

© rH 

O'! 

o 

CD 


00000©ffl00000000 

j-OOOOlDfflOOOOOOOO 

-OOOOOOOOOOOOOO 

iOHOrtHfOlDolDrlrlOOO 


HOfl-om 
O O O O O 

AAA AA 

O O ID O CD 
© O 00 O ID 
O o 5l O O 
O O ffl O H 

H o ffi O H 


(\J M H O O 

O O o o o 

AAA A A 

© rH O O O 
OTOOO 

© O O o 

rt ffl O O O 

rH <7 rH O O 


OOOOOOOOOOOOOO 


at 

o 

<r ~ 

!x 

3 at 
3 


o o o o 
O O O O 

* AAAA 

CL O O O O 

H O o o o 
<E o o o O 
H o o o o 
CD O O o o 


o o <-t «-t 
o o o o 
+ + i t 
© © © © 
O O <7 
o o o o 
O O <7 T 
O O CD CD 
O O CO CD 


rH rH rH O O O 

0 o o o o o 

1 i r + + t 
QCQ6C6 
<7 <7 <T O O O 
O O O O O O 

<7 <7 o o o 

CD CD CD O O O 
CD CD CO o O O 


CL 

£ 

o 

£ 

LJ 

~7 

5 

CL 


oooooooooooooo 


o 

u 

o o 
\z ° 
22 


rH rH «rt O rH rH O 

O O o o o o o 

£££££££ 

o o o a © o 

O O O o o o O 

OOOO^OO 
00 ffl CD O ffl CO O 

H N H 7 O 


rH 0J (\l rH ,rt O 

o o o o o o 

AAAAAA 

o © <H o O O 
O CD <7 O O o 

o © *7 o o o 
oo «-* © o o o 

rH rH «H H O 


oooooooooooooo 


>0000000000000 


LJ (DCDCDCD^CD^rrr^^^rCDCDCD 

at oooooooooooooo 

X ) * 

3 h- CD i 

<i<niNN 


OOOOOOOOOOOOOO 

)AA AAAAAAAAAAA A 

_NN<M<M©lX>0(MCJ©©fM©© 


r © in 

ID 

ID LI 

w-4 

66 cm 



rH 

f rl 

\r> co io 

LJ 


* 1 

* 

■ « 

■ 

• 


• 

* 

■ * 

H O 

O 

o o 

o 

o o 

o 

o 

o 

o 

o 

o o 

CJ 

CM 

CM CJ 

CD 

(D CD 

CM 

CJ 

CJ 

CJ 

CJ 

CJ CJ 

O 

O 

O O 

o 

O o 

o 

O 

o 

o 

o 

o o 

LJ 4* 

J at CM © 

£ 

AA 

+ 

© 

AA 

+ 

© 

+ 

© 

£ 

+ 

© 

A 

+ + 
© © 

<1 3 CL LO 

<0 

rH rH 

(0 

O CM 

ID 

ID 

CO 

io 

rH 

(M CJ 

|- © H CT) 

© 

CD CD 

o 

O N 

ID 

ID 

rH 

rH 

CD 

CJ CM 

o © <x © 

© 

<7 <T 

CD 

O CD 

CM 

CM 

iO 

CO 

<7 

rH H 

h LJ h 7 

<7 

ID ID 

ID 

ID CD 

sr 

<r 

CO 

CO 

ID 

ID ID 

a: cn rH 

rH 

rH rH 

rH 

rH rH 

CD 

CD 

rH 

rH 


W rH 

Q. 

• 


» 

• • 

• 

• 

• 


• 

• ■ 

O 

o 

o o 

O 

o o 

o 

o 

o 

O 

O 

o o 


gi 

O O 

I- o 

2g 


rHOOOOrHOrHrHrHrHrHO 

OOOOOOOOOOOOO 

AAAAAAAAAAA AA 

ooooooooNoooo 

000000©0ffl©000 

OOOOOOOQfflCMOOO 

OOOOOOOOCSJOOOO 

CJOGOOrlOHrlHrlHO 


OOOOOOOOOOOOOO 


OHOHHGHOH^rtHHH 

OOOOOOOOOOOOOO 

t AAA AAAAAAAAAAA 

boooooooooNMooo 

OoOOOOOOOOOCDOOO 

HOOOOOOOOONNOOO 

^oooooooooooooo 

OOHOH(MOhDHH(Mhh*H 


oooooooooooooo 

t rt AA AAAAAA AAAAAA 

XZCLOOOCMOvJCMOOOOOOCDCD 

OQ(-ooo0)fl)O(DCDC00!)C0oOO 
►-.J<TOO OCDCDCDC , Jf\JC , JC v 4CMVDkDCD 
Ul^l-ffiOCD^^^lOLDUDliDlDCDCDCD 
3 0)NNNCM<\ICM<MWfVIWNIDtDin 

oooooooooooooo 


OrHCVJCDT 

J^^rHDJCD^IDUlNCOCTt^l^rH^^ 

mCL^C\jCD’^-lDLDKCD0^Or-i(\jCD'7 
u H H H H rl 


oooooooooooooo 


o 

t i 

b o 
O o 
H O 


W O CD 

o o o 

AAA 

O O H 

CO o to 
10 ID O 
CM ID «H 
C\J rf ID 


HOP) 

0 o o 

1 + + 
O Cl Ci 
O O <7 
O o ID 
O O CD 
O O CD 
CM O ID 


O CD 

o o 

AA 

o o 
o o 
o o 
o o 
O CJ 


CM CJ O 

O o o 

AAA 

<7 (T\ o 
CM o o 
rt MO 
<tho 
00 ^ o 


O CD 

o o 

AA 

o o 
o o 
o o 
o o 
O CM 


*0000000000 


fcSI 

h6 

o 2 . 


N 

LO 

CD CO 
CO 


a 

(3 <x 
<r at 
at Cl 



at 

3 

_ 


X 2 


3 


\ s 

<1 


if) 

CD 3 

h 

+ 

LJ 

3 L. 

{- 


CD 

at 

<L 

3 

CD 

<-> 

O 

CD 


<T 

CL 


• <L CO <1 
\ \ ~i 
I H -j 
. tn « 
3 CL 
I CL O CD 

I X o 

i H . 


CD 

H 


Z 

3 

ID 

D 

z 

O 

n 


ID 

H 

CD 

ID 

<X 

CD 


at 

» 

CJ 

3 

ID ID 


h- 

(M O 



rH ID 



ID 


CO 




OJ CL 
X 


II 

> 

at 

LJ 


u 

LJ 

at 


to 

3 

at 

X 

h- 


if) 

U 

*: 

<r o 

at L. 
CD CO 
h- 
_J 

<1 © 
K U 
O _l 
H -J 

<r 

H 

cn 

z 

o 


II 

LJ 

o 


H 

<L 


CO 

H 


ID 

CM 


ID 

CM 


u 

LJ 

CO 

\ 


at 

at 

at 

U 


3 


3 






o 

© 

H 

if) 

U 

CD 

X 

CO 

z 

CO 

z 

3 

3 

H 

1“ 


• 

-J F; 

3 

hi 

|- 

U 

3 





3 

3 

L 

3 

© 

o 

w CD 

3 

at 

P 

etc 


\— 

© 

H 

m 

N 

N 

<x 

5 



3 Z 

CL 


CL 

H 

L) 

at 

U 

a: 

N 

N 

X 

I 

o 


3 at 

n 

r 

3 

r 

C 

3 

3 

3 

3 

O 

O 

ID 

CD 

3 

If 

O X 


5 

Cl 

o 

LJ 

o 

H 

© 

H 

Z 

Z 




X 

3 K 

3 

o 

CD 

a 

X 










o 

3 

<x 















at (- 

H 

CJ 

CD 

<7 

ID 

cr* 

N 

CO 

© 

o 

CM 

CD 

^7 

ID 

r 

« 3 

O 









rH 

rH 

rH 

rH 

rH 


-X Z 



Sfc 

a : 3 
© a: 
X 


CD 

if) 

ctf 

C 3 

CD 

a 

E 

ca 

CO 

u 

o 

L~l 

CD 


lx. 


o 

Cl 

4 -> 

3 

o 


>> 

u 

CD 

c 

•-H 

x : 

u 

cti 

E 

o 

X 

Lj 

3 

{—• 


iT) 

\o 

CD 

u 

3 

60 

Uh 


27 



The brake specific fuel consumption stays relatively constant down to about 
657. of full power and then starts to increase. This can be attributed to the 
influence of the power recovery from the exhaust turbine which tends to keep 
the BSFC low even at part load. 


6.2 Case Number Two: Opposed Piston 3 Cylinder Diesel Engine 

This case is similar to the combined-cycle engine described in case 
number 1. The purpose of the engine and its design requirements are the same. 
The diesel input data file, shown in Figure 6.6, is the same as for the V-6 
engine except for the volume mode indicator that denotes an opposed piston 
configuration. Also, the exhaust port width is specified instead of the 
exhaust valve area profile. The turbomachinery input file shown in Figure 6.7 
is the same as the earlier case except that the fuel injected per cylinder is 
doubled to reflect the increased air flow provided per cylinder. 

The design point output files are shown in Figures 6.8 and 6.9, 
respectively. In general, the differences between the two engines are small. 
There are almost no differences in the turbomachinery. The most significant 
differences in the diesel engine are the peak pressure and the equivalence 
ratio. The difference in peak pressure probably results from the difference 
in effective compression ratio that occurs when the exhaust port of the 
opposed piston engine closes 21 degrees later than the exhaust valve in the 
V-6 engine. The difference in equivalence ratio is probably due to 
differences in scavenging that also relate to changes in valve/port timing. 

No attempt was made to optimize the port timing. Performance improvements 
beyond that shown may be possible. 
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TW 0 - ST R 0 K E , OF P 0 S E D 

3.100 2.340 7.168 

1 , 5 345 . 65 . 0 . 0 

0.55 1460. 1460. 

0 100. 260. 0 1 
0.5542 0.8 

0 . 5542 0 . 8 


PISTON 3- 

3.1712 
3 . 5 
1460 . 


CYL 


Figure 6.6 


Diesel Input Data File for 


ENG I NE 


Sample Case 2 
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TEST CASE: COMBINED CYCLE ENGINE 

SD I CEC=0 ,NCODE=l , LONG=T , NMODES=l , SEND 
SD MODE=l , 

KONFI G< 1 , 1)=1 ,1 ,0 ,2,0, SPEC d , 1 ) =7 . 3 , 4*0 , 1 , 5*0 ,0,0, 

KONFI Gd ,2) =4, 2,0 ,3,0 ,SPEC(1 ,2) =1 .8,0,1,1001,1 ,1002,1 ,1003,1 , 
2. ,0 , .35,1 .05,1 , 

KONFI G( 1 , 3) =7,3, 0 , 12,4 , SPECd , 3) =0 . 455 ,0 ,0 , 

KONFI G( 1 ,4) =4, 4, 0,5,0, SPECd , 4) =1 . 3 , 0 , 1 , 1004 , .5,1005,1,1006,1, 
0,0, .73,3.3132,1 , 

KONFI Gd , 5) =6 , 12 , 5 , 13 , 6 , SPECd » 5) = . 02 , . 02 , 250 . , 0 . 4 , 1 , 

KONFI G(1 ,6) =2, 6,0 ,7,0 , SPEC( 1 , 6) =0 .1 ,0 ,0 .00023,6122. ,1,13300. , 
FARRAYd ,6) =2,1 

KONFI Gd ,7) =5,7 , 0 , 3 , 0 , SPEC( 1 , 7) =4 . 3 , 0 , 1 , 1007 , 1 , 1003 , 1 ,1 , 

0,1, .845,5000,1, 

KONFI G(l,3)=2, 3, 0,3,0, SPEC (1,3) =0,0, 0,0, .33,13300 . , 

FARRAYd ,8) =2,1 

KONFI G(1 ,3) -3,9,0,10, 0,SPEC( 1,9) =1.8, 0,1 ,1009,1 ,1010,1 ,1 , 

0,1, .343,5000,1, 

KONFI G(1 ,10) =9,10 ,0 ,11 ,0 , SPECd ,10)=0 . ,1 ,0 . ,0 . , .3350,1. ,0 . , 

0 . , 1 . , 

KONFI Gd ,12) =3, 13, 0,14,0, SPEC< 1,12) =0,1, 0,0. ,0.3850,1 . ,1 . , 

0 . , 1 . , 

KONFI G( 1,13) =11, 2, 4, 7,0, SPEC ( 1 , 13) =3*1 . 

KONFI G(1 ,14) =11, 3, 15, 0,0 ,SPEC(1 ,14) =3*1 . 

KONFI G( 1,15) =10, 0,0,0 , 0 , SPEC( 1 , 15) =-200 , 0 , 0 , 

KONFI GC 1 ,16)=12,SPCNTL(1 . 16) =1 , 9 , 100 , 3 , 10 , 0 , 0 , 

KONF I Gd ,17) =12, 3PCNTL( 1 ,17) =1,7, 100 ,3, 3, 0,0, 

KONFI G( 1 ,13) =12, SPCNTLC1 ,13)=1 ,4,100 ,3,7.0 ,0 , 

KONF IG<1 ,20)=12,SPCNTLd ,20) =1 ,2,100 ,3,4,0 ,0 , 

KONFI G(1 ,21)=12,SPCNTL(1 , 21 ) =1 ,1 ,100 ,3,2,0 ,0 , 

KONFI Gd , 22) =12 , SPCNTLC 1 , 22) =1 , 13 , 200 , 3 , 13 , 0 , 0 , 

KONFI G<1 , 23) =12 , SPCNTLC 1 , 23) =1 , 15 , 200 , 3 , 14 , 0 , 0 , 

KONFI G(1 , 24) =12 , SPCNTLC 1 , 24) =3 , 5 , 200 , 3 , 5 , 0 , 1 , 

KONFI G(1 ,26) =12, SPCNTL (1,26) =10, 2, 20 0,5, 2, 4. 0,0, 

SEND 

SD SPEC( 3 , 16) =1 , SPEC< 3 , 17) =1 , SPEC< 3 , 13) =1 , 

SPEC( 9 , 20 ) =1 , SPECC 3,21) =1 , SPECC 9 , 22) =1 , SPEC( 3 , 23) =1 , 

SPECC 3 , 26) =1 , SEND 
SD SPECC 3, 6) =0.00024 SEND 
SD 3PEC<3,6)=0 .00020 SEND 
SD SPEC<3,6)=0. 00013 SEND 
SD SPEC<3,6)=0 .00016 SEND 
SD SPEC( 3, 6) =0.00014 SEND 
SD SPEC( 3, 6) =0.00012 SEND 
SD SPEC(3,6)=0.00010 SEND 
SD SPEC(3,6)=0 .00009 SEND 
SD SPEC< 3, 6) =0.00008 SEND 
SD SPEC(3,6)=0 .00007 SEND 
SD SPEC(3,6)=0 .00006 SEND 


Figure 6.7 Turbomachinery Input Data File for Sample Case 2 
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TWO-STROKE, OPPOSED PISTON 3 

BORE 

STROKE 

CONNECTING ROD 
GEOMETRIC COMPRESSION RATIO 
EFFECTIVE COMPRESSION RATIO 
NUMBER OF CYLINDERS 
SWEPT VOLUME (TOTAL) 

ENGINE SPEED 
MEAN PISTON SPEED 
SUPPLY AIR PRESSURE 
SUPPLY AIR TEMPERATURE 
EXHAUST PRESSURE 


-CYL 


7 

9 

6 

3 

133 

6122 

3000 

150 

885 

129 


ENGINE 
.100 INCHES 
.940 INCHES 
.163 INCHES 
.17 
.21 

.141 CUBIC INCHES 
. RPM 
. FT/M IN 
.00 PSIA 
.86 DEG R 
.66 PSIA 


PORT TIMING: 


INTAKE OPEN 125. 
INTAKE CLOSE 235. 
EXHAUST OPEN 100. 
EXHAUST CLOSE 260 , 


FUEL FLOW RATE 
AIR FLOW RATE 


0.0833 LBM/SEC 
2.4390 LBM/SEC 


FUEL INJECTED/CYCLE 0.000290 
AIR INDUCTED/CYCLE 0.007963 


IMEP 

BMEP 

FMEP 


421 .4 PSI 

393.4 PSI 
23.0 PSI 


I HP 
BHP 
FHP 


867.3 HP 
810.1 HP 
57.7 HP 


I SFC 
BSFC 


0.3633 LBM/HP-HR 
0.3945 LBM/HP-HR 


INDICATED THERMAL EFFICIENCY 
BRAKE THERMAL EFFICIENCY 


0 . 3795 
0.3542 


MAXIMUM PRESSURE 
MAXIMUM TEMPERATURE 
MAXIMUM RATE OF PRESS RISE 


2679.4 PSIA 

4662.5 DEG R 
62.01 PSI/DEG 


AT 374.5 DEG 
AT 337.4 DEG 
AT 359.4 DEG 


PURITY 0.3287 
SCAVENGING EFFICIENCY 0.S058 
TRAPPING EFFICIENCY 0.6096 
CHARGING EFFICIENCY 0 .4138 
DELIVERY RATIO 0.6733 
RESIDUAL FRACTION 0.1973 


EQUIVALENCE RATIO DURING COMPRESSION 0.1620 
EQUIVALENCE RATIO AT EVO 0.3530 
EQUIVALENCE RATIO BASED ON FUEL AND 

AIR FLOW RATES (EXHAUST) 0.5233 


% FUEL 

ENERGY 

TO 

HEAT 

LOSS 

4.70 

% 

% FUEL 

ENERGY 

TO 

WORK 


37.95 

% 

% FUEL 

ENERGY 

TO 

EXHAU 

ST 

57 . 35 

% 

EXHAUST 

TEMPERATURE 

2187 

.3 DEG 

R 


Figure 6.8 Diesel Output File for Sample Case 2 


LBM/ CYCLE 
LBM/CYCLE 
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CA PCYL TCYL PHICYL MCYL MDOTIN MDOTOUT BALANCE 

100.00 405.4 3378.9 0 . 8580 0.006350 0.000000 0.000000 0.00000 

110.00 325.2 3221.2 0.8530 0.005904 0.000000 0.000081 0.00000 

120.00 230.3 2990.1 0.8580 0.004890 0.000000 0.000113 0.00000 

125.00 188.2 2882.4 0.8580 0.004325 0.000000 0.000110 -0.00003 

130.00 153.4 2737.5 0 . S580 0.002139 -0.000005 0.000085 -0.00003 

140.00 129.7 2439.2 0.7413 0.002454 0.000054 0.000008 -0.00003 

150.00 130.1 2138.4 0.5741 0.003134 0.000080 0.000018 -0.00004 

160.00 130.8 1883.8 0.4417 0.004038 0.000099 0.000027 -0.00001 

170.00 131.0 1694.0 0.3485 0.005084 0.000109 0.000035 0.00000 

180.00 131.3 1558.5 0.2846 0.006199 0.000112 0.000039 0.00009 

185.45 131.5 1502.3 0.2587 0.006808 0.000111 0.000041 0.00011 

190.00 132.6 1465.4 0.2406 0.007010 0.000105 0.000069 0.00026 

190.00 132.5 1465.5 0.2407 0.007010 0.000105 0.000069 0.00010 

200.00 133.3 1395.9 0.2090 0.007278 0.000093 0.000075 0.00010 

210.00 133.7 1344.6 0.1862 0.007371 0.000074 0.000073 0.00009 

220.00 134.0 1310.2 0.1710 0.007236 0.000049 0.000067 0.00009 

230.00 134.1 1292.2 0.1630 0.007014 0.000018 0.000054 0.00004 

240.00 134.6 1291.4 0.1620 0.006590 0.000000 0.000041 0.00004 

250.00 139.5 1304.1 0.1620 0.006230 0.000000 0.000030 0.00005 

260.00 153.7 1337.8 0.1620 0.006055 0.000000 0.000000 0.00005 

270.00 180.6 1395.1 0.1620 0.006055 0.000000 0.000000 0.00005 

280.00 218.2 1464.8 0.1620 0.006055 0.000000 0.000000 0.00005 

290.00 272.0 1549.4 0.1620 0.006055 0.000000 0.000000 0.00005 

300.00 350.4 1651.3 0.1620 0.006055 0.000000 0.000000 0.00005 

310.00 467.0 1773.2 0.1620 0.006055 0.000000 0.000000 0.00005 

320.00 641.5 1916.2 0.1620 0.006055 0.000000 0.000000 0.00004 

330.00 896.2 2077.0 0.1620 0.006055 0.000000 0.000000 0.00004 

340.00 1233.0 2240.4 0.1620 0.006055 0.000000 0.000000 0.00004 

350.00 1602.0 2413.7 0.1698 0.006058 0.000000 0.000000 0.00004 

360.00 2152.1 2997.4 0.2748 0.006102 0.000000 0.000000 0.00007 

370.00 2633.6 3890.5 0.4889 0.006191 0.000000 0.000000 0.00012 

380.00 2569.9 4518.5 0.6982 0.006278 0.000000 0.000000 0.00010 

390.00 2091.8 4658.4 0.8140 0,006327 0.000000 0.000000 0.00013 

400.00 1567.7 44S9.7 0.8510 0.006342 0.000000 0.000000 0.00012 

410.00 1164.8 4237.9 0.8576 0.006345 0.000000 0.000000 0.00012 

420.00 886.5 4002.8 0 . 85S3 0.006345 0.000000 0.000000 0.00012 

430.00 696.6 3802.9 0.8584 0.006345 0.000000 0.000000 0.00011 

440.00 565.1 3635.7 0.85S4 0.006345 0.000000 0.000000 0.00011 

450.00 472.2 3496.6 0.8584 0.006345 0.000000 0.000000 0.00011 


Figure 6.8 - - continued 
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7 Conclusions 


The computer program described in this report can simulate both design 
point and off-design operation of a combined-cycle engine. It retains the 
flexibility of the two original programs on which it is based. Widely varying 
engine configurations can be evaluated for their power, performance and 
efficiency as well as the influence of altitude and air speed. Although the 
program was developed to simulate aircraft engines, it can be used with equal 
success for stationary and automotive applications as well. 

The program has been tested for a wide variety of engine configurations. 

A student project is currently underway at Iowa State University to use the 
program to simulate potential engines for a high-altitude, long-endurance 
aircraft such as the Boeing Condor [8]. This is the type of application for 
which the program was intended and it should provide useful information to 
allow comparisons between engine configurations as well as for optimizing a 
given engine design. 
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9.2 Appendix B: Second Law Analysis 

Second law analysis is a technique for analyzing thermodynamic systems 
that goes beyond traditional energy balances based on the first law of 
thermodynamics. Second law analysis uses the second law of thermodynamics to 
determine the fraction of the energy stored or transported in a process that 
could be converted to useful work. This type of analysis uses a property 
called "availability." 

Availability is the maximum theoretical work output that could be 
developed in bringing a system into equilibrium with a reference environment. 
In general, this involves all work that could be obtained by exploiting 
differences in velocity, potential energy, temperature, pressure, and 
chemical potential between the system and the environment. When thermal, 
mechanical, and chemical equilibria are attained and the system is at rest and 
zero elevation relative to the reference environment, no further potential 
exists for developing work. In this condition, the system and environment 
combination is said to be at the dead state. The value of the availability 
depends on both the state of the system and the state of the environment; it 
is a measure of the system's thermodynamic potential relative to the chosen 
environment L 9 J . 

Availability is not a conserved quantity. It can be destroyed by 
irreversibilities within the system such as friction, heat transfer, mixing, 
turbulence, etc. However, it is possible to combine these terms into a single 
loss term and write an availability balance for a system that is analogous to 
the energy equation. 

dA l „ dV\ m' r . 

-WO-rvTKM- (»"-*■-*) (9 u 

1 2 3 4 5 

The numbered terms in this equation can be interpreted as availability 
quantities as follows: 

Term 1: the rate of change of availability in the control volume 

Term 2: the rate of availability transfer associated with heat transfer 
from the gases to the wall 

Term 3: net rate of availability transfer to the piston associated with 

work 

Term 4: rate of availability input to the control volume with the fuel, 

including the availability input associated with flow work 

Term 5: rate of availability destruction due to irreversibi 1 it les 

within the control volume 

Second law analysis consists of applying the availability balance to the 
components of a system to determine whether the processes in that component 
are destroying availability. This allows the individual processes within a 
system that responsible for inefficiency to be identified. 
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Computer Program ACV 


A computer program, called ACV, has been written that performs a second 
law analysis of a combined cycle engine. The program is interactive, and the 
user enters information such as temperature, pressure, mass flow rate, etc., 
from the output files of the cycle simulation program, CCE. The output is 
given in terms of the availability transferred due to heat transfer, work, 
and mass flows as well as the availability destroyed by irreversibilities. 

A sample output from the program is shown in Figure 9.1. This case 
corresponds to the V-6 engine sample case discussed in Section 6.1. The 
results can be summarized as: 


Availability input with air 

1.68 

Availability input with fuel 

1796.26 

Total availability added: 

1797.94 Btu/sec 

Work 

711.54 

Heat 

87.89 

Irreversibi lites 

637.34 

Availability exiting with exhaust 

342.74 

Availability exiting with air 

18.63 

Total availability accounted for: 

1798.14 Btu/sec 


These results show that of the availability added to the combined cycle 
engine, 711.54/1797.94 = 39.67. is converted to work. This compares to a first 
law brake thermal efficiency of 44.07. based on the lower heating value of the 
fuel. The analysis also shows that the most significant loss is the 
irreversibility of the combustion process. 27.67. of the availability is 
destroyed by the irreversible combustion process, 3.64 is destroyed in the 
heat exchanger and 19.17 exits with the exhaust gases. 
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Figure 9.1 Sample Output from Second Law Analysis Program ACV 



9.3 Appendix C: Use of Centrifugal Compressor Data in the Combined- 
Cycle Engine Simulation 


The computer program NNEPEQ is accompanied by a set of default compressor 
and turbine maps that are scaled from a 100 lb/sec axial-flow compressor and 
turbine. Since a 1000 hp combined-cycle engine such as that discussed in 
Section 6.1 is likely to use a 2 Ib/sec centrifugal compressor, the effect of 
using actual data for a compressor of this type was investigated. Skoch and 
Moore [10] collected experimental data on the performance of two 10 lb/sec 
centrifugal compressors. The blade thickness of one of the compressors was 
increased to simulate a scaled-up 2 lb/sec centrifugal compressor. While the 
performance data for the simulated 2 lb/sec compressor would have been ideal, 
only data for the thin-blade compressor were available. These data were used 
to create the compressor map shown in Figure 9.2. When this map was used for 
the high pressure compressor in a sample case of a combined-cycle engine, the 
results shown in Table 9.1 were obtained. 


Table 9.1 Comparison of Default Compressor Map and 10 Ib/sec Centrifugal 
Compressor Map 


rbine Inlet 
emperature 



Default 

Compressor 


Centrifugal 

Compressor 



(deg R) 










HP 

Surge Mar 

Effic 

PR 

HP 

Surge Mar 

Effic 

PR 

2300 

454 

30.2 

78.9 

10.1 

455 

49.1 

79.0 

10.1 

2200 

417 

34.4 

80.4 

9.6 

394 

51.9 

78.0 

9.1 

2100 

377 

37.1 

81.6 

8.9 

344 

54.9 

77.9 

8.3 

2000 

335 

38.3 

82.5 

8.3 

284 

57.6 

76.7 

7.2 

1900 

286 

37.2 

82.6 

7.4 

232 

59.7 

. 76.0 

6.3 

1800 

224 

33.9 

80.6 

6.2 

187 

61.9 

75.8 

5.4 

1700 

170 

31.6 

79.0 

5.2 

134 

63.3 

74.6 

4.4 

1600 

123 

25.6 

78.0 

4.2 

88 

63.2 

74.0 

3.5 

1500 

78 

9.1 

77.0 

3.3 

56 

60.6 

75.7 

2.8 


As is clear from the table, the centrifugal compressor maintained a wider 
surge margin and higher efficiency than the default compressor. The pressure 
rise through the centrifugal compressor dropped off somewhat faster as the 
turbine inlet temperature was decreased. Although not shown, a second series 
of cases was also run with the "R" value of the centrifugal compressor 
adjusted so that the initial surge margin of the default compressor and the 
centrifugal compressor matched. In this case, the relationship between the 
turbine inlet temperature and the compressor pressure ratio and efficiency was 
altered somewhat but the performance was still below the default compressor. 
This study makes it clear that while the default maps can be used for 
comparisons between configurations, actual maps are required to obtain 
accurate performance predictions. The default maps were used for the sample 
cases discussed in Section 6.0. 
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Figure 9.2 Centrifugal Compressor Map 



9. A Appendix D: Subroutine for Simulating Opposed Piston Engines with 
Crankshafts Out-of-Phase 


A FORTRAN subroutine was developed for use with opposed piston engines 
with out-of-phase crankshafts. The subroutine has been tested but was not 
incorporated into the CCE simulation program. Since this type of engine is 
very specialized, the added difficulty to the program user of having to input 
additional engine specifications was considered undesirable. The subroutine 
could be added with only minimal program modifications. The volume subroutine 
listing, along with a simple test program, is provided in Figure 9.3. 
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ono o non o on nnnon 


THIS ISA PROGRAM TO TEST THE VOLUME SUBROUTINE FOR AN OPPOSED 
PISTON ENGINE WITH CRANKSHAFTS OUT-OF-PHASE 

COMMON/ENGINE/ CONROD , BORE , STROKE , CR , VCL , MODE , NCYCLE , LTV PE 
COMMON/OPPSED/ PSI 
SPECIFY VALUES OF INPUT VARIABLES 
C0NR0D=9. 

B 0 R E = 4 . 5 
STROKE-4 . 5 
CR=1 6 . 

MODE-1 
PSI=10 . 

CALCULATE CLEARANCE VOLUME (REQUIRED BEFORE FIRST CALL TO 
VOLUME SUBROUTINE) 

CALL CLEAR 

WRITEC5,*) 'VCL = ' ,VCL 
DO CO 1=1,21 
CA=-1 0 . +FLOAT ( I -1 ) 

CALL V0L2( CA , V , DVCA) 

WR I TE ( 6 , * ) CA , V , DVCA 
20 CONTINUE 
STOP 
END 

SUBROUTINE V0L2( CA , V , DVCA) 

COMMON/ENGINE/ CONROD , BORE , STROKE , CR , D , VCL , MODE , NCYCLE , LTYPE 
COMMON/OPPSED/ PSI 

CALCULATE CYLINDER VOLUME AND VOLUME DERIVATIVE 

DUMMY 1 =SQRT ( CONROD** 2- ( STROKE/ 2 .*SIN< CA/57 . 29578) >**2) 
DUMMY2=SQRT(C0NR0D**2-< STROKE/2. *S I N( ( CA+PSI )/57 . 29578) )**2) 
XL= STROKE/ 2 . * ( 1 . -COS ( CA/57 . 29578) )+CONROD-DUMMYl 
XR= STROKE/ 2 . ★ ( 1 . -COS< (CA+PSI )/57 . 29578) )+ CONROD -DUMMY 2 
V=VCL+0 . 7853982*B0RE*B0RE*(XL+XR) 

DXLDT=STR0KE/2 .*SIN< CA/57 . 29578)+ 

1 ( 1 .+STR0KE/2 . *COS < CA/57 . 29578) /DUMMY1 ) 

DXRDT =STR0KE/2 .★SI N( ( CA+PS I )/57 . 29578 )★ 

1 (1 ,+ STROKE/ 2 , ★COS ( ( CA+PSI )/57 . 29578 ) /DUMMY2 ) 

DVCA=0 . 785393 2^ BOR E+ BOR E+( DXLDT+ DXRDT ) /57 , 29578 

RETURN 

END 

SUBROUTINE CLEAR 

COMMON/ENGINE/ CONROD , BORE , STROKE , CR , VCL , MODE , NCYCLE , LTYPE 
COMMON/OPPSED/ PSI 

CALCULATE CYLINDR CLEARANCE VOLUME 


CA=-PSI/2. 

DUMMYl=SQRT(CONR0D**2-( STROKE/2 .★SIN (CA/57. 29578) )**2) 
DUMMY2=SQRT(C0NR0D**2-< STROKE/2 .★SI N( < CA+PSI )/57. 29578) >**2) 
X L = ST R 0 K E/ 2 . ★ ( 1 . -COS ( CA/57 . 29578) ) + CONROD- DUMMY 1 
XR=STR0KE/2 .★( 1 . -COS( ( CA+PSI )/57 . 29578) ) +C0NR0D-DUMMY2 
VM I N = 0 . 7853982*B0RE+B0RE+(XL+XR) 

CA=CA+1 80 . 

DUMMY1 =SQRT(C0NR0D**2-(STR0KE/2.*SIN( CA/57. 29578) )**2) 
DUMMY2=SQRT ( C0NR0D**2-( STROKE/2 .★SIN ( ( CA+PSI )/57. 29578) )**2) 
XL= STROKE/2 .★( 1 . -COS( CA/57 . 29578) )+CONROD-DUMMYl 
XR=STR0KE/2 . ★ ( 1 . -COS ( (CA+PSI )/57 . 29578 ) ) +C0NR0D-DUMMY2 

VMAX=0 . 7853982*B0RE*B0RE*(XL+XR) 

VCL= ( VMAX-CR+VMI N)/( CR-1 . ) 

RETURN 

END 


Figure 9.3 Cylinder Volume Subroutine for an Opposed Piston Engine 
with Out~of “Phase Crankshafts 
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